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Approximating Ropelength by Energy Functions 

John M. Sullivan 



Abstract. The ropelength of a knot is the quotient of its length by its thick- 
ness. We consider a family of energy functions for knots, depending on a 
power p, which approach ropelength as p increases. We describe a numerically 
computed trefoil knot which seems to be a local minimum for ropelength; there 
are nearby critical points for , which are evidently local minima for large 
enough p. 



1. Thickness and ropelength 

We measure the ropelength of a knot as the quotient of its length by its thick- 
ness. The thickness is the radius of the largest embedded normal tube (called the 
thick tube) around the knot. Ropelength is a mathematical model of how much 
physical rope it would take to tie the knot. In |CKSOl| we showed that in any knot 
or link type there is a ropelength minimizer, and that minimizers are necessarily 
C^'^ curves. But still, explicit examples of tight (ropelength minimizing) links are 
known only in very special cases. 



One way to define the thickness of a link (following [GM99]) is to consider 
circles through three points on the link. For any three distinct points x, y, z in R'^, 
we let r{x,y,z) be the radius of the (unique) circle through these points (setting 
r = (X) if the points are collinear). Also, if is a unit vector at x, we let r{Tx, y) 
be the radius of the circle through y tangent to &t x. Note that 



■r{Tx-,y) = 



2^\\x-y\\^-{n-{x-y)y 



Now let L be a link in K^, that is, a disjoint union of simple closed curves. We 
define the thickness t{L) of L as 

t{L) := inf r{x,y,z). 
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This quantity vanishes when L fails to be C^'^, so from now on we wiU restrict 
attention to C^'^ hnks. 

The infimum in the definition of r is always achieved when (at least) two of the 
three points approach each other |GM99, [CKSOl : we have 



t{L) = inf r{Tr,L,y), 

y=^x 

where T^L denotes the tangent vector to the curve L at x. We can equiva- 



lently |CKS01] define the thickness in terms of the medial axis (see [A.CK01|) 
of L. The medial axis is the set of points in space which fail to have a unique 
nearest point in L. Then the thickness t{L) is exactly Federer's reach [ Fed59| 
of L, the distance from L to its medial axis. 

Gonzalez and Haddocks [ |GM99| considered the quantity 

T^^{L) inf r{x,y,z) = inf r(TyL, x). 

v,^ V 

which they called the "global radius of curvature" but which could well be termed 
the local thickness of the link L at the point x E L. However, it seems more natural 
to define local thickness in terms of higher-order contact at x: 

T.j,{L) := inf r{TxL,y). 

Either definition is a local thickness in the sense that 

r(L)=infT,(L) = infrG^(L). 

X X 

However, Tx{L) has a nice geometric interpretation as the radius of the largest 
sphere which, staying tangent to L at x, can be rotated completely around L 
without touching L at any other point. 

Given a circle C of radius r, consider the set of points in space whose distance 
to C is less than r. Following Dan Asimov, we will call this a hialy (or more 
precisely an open solid bialy), with neck at the center of C and axis normal to C. 
A bialy is an open solid torus of revolution whose major and minor radii are equal. 
Given a point a; on L, the bialys with neck at x and axis along T^L are nested, and 
it is clear that Tx{L) is radius of the largest one which avoids L. 

We now check that our definition of local thickness also agrees with the notion 
of local feature size from the theory of medial axes. Let r{T,j;L,TyL) denote the 
radius of the sphere tangent to both T^L and TyL (the smallest such sphere if these 
vectors are cocircular). Then it is clear from our discussion above that 

Tj;{L) = inf r{TxL,TyL), 

V 

since the interiors of all spheres of radius r tangent at x fill out the bialy at x. 

A formal definition of medial axis is the following. Suppose L is a compact 
subset of R'^. Then for any point p G M"^, we can consider the distance to p as a 
function dp : L — > R. It achieves its minimum d{p, L) on some set min{dp) C L. 
Then the medial axis is M{L) := {p E M.^ : #min((ip) > 1}, the set of p for which 
this minimum is achieved at more than one point. At a point p G A/, the local 
feature size is simply the distance d{p^ L) to L. At a point x E L, the local feature 
size is 

lfs(a;) u\i{d{p,L) : p E M,x <E m\n{dp)}. 
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(We note that this definition, intended for use when L is a curve or surface, is only 
one of many inequivalent definitions for local feature size that have been given in 
the computer science literature.) 

In terms of these definitions, we are now ready to state our first proposition: 

Proposition 1.1. The local thickness Tx{L) of a link L is also the local feature 
size of L at X. 

Proof. The medial axis M{L) is exactly the set of centers p of spheres which 
are tangent to L at two or more points, and whose interiors avoid L. The local 
feature size at p e A/ (or at the points of tangency along L) is the radius of the 
corresponding sphere. But the local thickness at x is also the radius of the smallest 
bitangent sphere tangent at x. □ 



2. Smooth approximations to ropelength 

We can approximate the infimum in the definition of ropelength with ener- 
gies. For a link L of length 1, we define 



R^iL) / riTxL,yyPds, i?f(L) f / R^^L) ds 



i/p 



We extend this energy to links of arbitrary length so that it is scale-invariant. Aside 
from this rescaling, the energy Rp{L) is the same as the energy Up,2{L) suggested 



in |GM99| 



Then clearly Rp{L) is bounded above by the ropelength R{L), and in fact 

lim i?P(L) = R{L). 

For p > 1, it also follows from Holder's inequality that Rp{L) > R^{L). The 



examples in |CKS01] show that tight links (ropelength minimizers) need not be 
smoother than C^'^. However, since the energies RP smooth out the hard-shell 
potential involved in a thickness constraint, we expect that minimizers for the 
energies Rf should be smooth. 



Diao, Ernst and Janse van Rensburg |DEJvR98| have listed several properties 



that one might desire for knot-energy functionals. For any p > 1, we conjecture 
that the energy Rp is "basic" in the sense that the absolute minimum is achieved 
uniquely by the round circle (for which R^ = 27r). It would suffice to prove this 
conjecture for the special case p = 1. 

For p < 2, the quantity Rp is not very useful as a knot energy, because it does 
not have an infinite barrier to self-crossings. It is straightforward to check that the 
energy of two skew lines which approach a right-angled crossing becomes infinite if 
and only if p > 2. 

The following conjecture would imply that, for p > 2, Rp is "charge" and 
"tight" , meaning that it approaches infinity for sequences of links which approach 
a curve with a self-intersection, or in which a knotted arc shrinks to zero size. 

Conjecture 2.1. For p > 2, the energy Rp is bounded below by some mono- 
tonic function of ropelength (so RP approaches infinity for any sequence of links 
with fixed length and thickness approaching zero). 

We will prove part of this conjecture, enough to conclude that RP is "charge" . 
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Proposition 2.2. If Lq is the smooth limit of smooth links Ln, and Lq has a 
self-crossing, then for p > 2 the energies RP{Ln) approach infinity. 

Proof. Let K/2 be an upper bound for the curvature of Lq, so that K bounds 
the curvature of Ln for large n. There is also a uniform bound T on the third 
derivatives. The thickness of L„ approaches 0. 

Let L be any of the L„ with n large enough that its thickness r satisfies 

r < ^, r « 1/T. 

Since r < 1 / K , the thickness is achieved by a doubly critical pair of points x, y, with 
L perpendicular to the segment xy at both ends. Picking an appropriate coordinate 
system, we can assume that x — (0,0, t) and y — (0,0, — r), with T^L = (1,0,0) 
and TyL = (cos </>, sin </>, 0) for some angle 0. We can then expand L in Taylor series 
near x and y as: 

a{s) = (s,/cis^r + fc2S^) + 0(/), 

(3{t) = {tcos(t) + ht'^ sm(t),tsm(j) - kst'^ cos(l),-T - k4t'^) + 0{t^), 

where the ki are components of the curvature vector, so < and the constants 
in the first omitted terms depend only on T. The tangent vector at x = a{s) is 
(1, 2fcis, 2fc2s) + O(s^). Thus for y = f3{t) we can compute 

r{T,L, y) = r + ^^(^L _ ^2) + t^{^{2 ~ sin^ + f ) - + fcs) + ©(s^, i^). 

This gives 

riT^L, y) < r (1 + C(|s| + \t\)^ + Ois^ t')) , 

where C = l/(4r^) + K/t. Since Kt < 1/4, we have C < 1/{2t'^). Now pick some 
small length a = t'^, and consider the integral defining Rp{L) taken over just the 
ranges < s,t < a. We get 



a pa 



RP{LY > / / r^Pdsdt 
Jq Jo 



T 



-P jj{l-pC{s + tf + 0{s^,t^)) dsdt 
-P {a'~pC{la') + 0{a')) 



> T-P (r29 - pr^Vr^ + C>(r^«)) 

Now if g > 1, the first term dominates the others. (We should note that a more 
careful analysis would confirm that the constants in the term of order r^'^ are still 
uniformly bounded in terms of T.) For 1 < g < p/2, this expression diverges to 
infinity as r — > 0, as desired. □ 

For p > 2, we expect that Rp is also "strong", in the sense of allowing only a 
finite number of link types under any given energy level, but this seems hard to 
prove. 

3. Numerical simulations 



We have earlier reported |KS97] on numerical simulations, conducted us- 



ing Brakke's Evolver |Bra92], of knots minimizing the Mobius- invariant energy 



of I FHW94 1 . We have recently used the Evolver for numerical simulations of tight 



(ropclcngth minimizing) links, using a discretization, for polygonal approximations. 



APPROXIMATING ROPELENGTH BY ENERGY FUNCTIONS 



5 




Figure 1. A pair of orthogonal projections (printed for cross-eyed 
stereo viewing) of the trefoil knot 7p for p = 4096. This is a numerically 
computed local minimum for the energy R^, and it is presumably close 
to a local minimum for ropelength. Its ropelength is about 37.5, as 
opposed to 32.7 for the minimizing trefoil. 



of the energies Rp introduced here. These simulations will be described in more 
detail in a forthcoming paper. 

As an example of the results, we describe here a trefoil knot, shown in Fig- 
ure ^ which seems to be a local mininum for ropelength. This trefoil is embedded 
in space as a (3, 2)-torus knot, as opposed to the (2, 3)-torus knot which gives 
the (presumed) global minimum. Its symmetry group in M'^ is 2 2 2 in the Con- 
way/Thurston notation for 2-orbifolds: that is, it is generated by the rotations of 
order 2 around three orthogonal axes. 

Numerically we have computed a family of very similar curves 7p with this 
symmetry which are critical points for the energies R^. Presumably, they are global 
minima if the symmetry is fixed. But without enforcing the symmetry, they seem 
to be unstable saddle points for Rp when p < 512 but stable local minima when 
p > 1024. (We have only tested powers p which are powers of 2.) In Figure || we 
plot the curvature and torsion of 740961 as functions of arclength along the curve. 

Presumably, the limit of these curves as p — > cx) is a curve 7 which is a local 
minimum for ropelength. This limit curve is interesting in that it seems to include 
two straight segments. The thick tube around these arcs does not contact that 
around any other part of 7. Presumably, at the ends of these segments, where 
the tubes do contact, the curvature of 7 jumps discontinously to some positive 
value (close to 1 if we normalize to t = 1). Thus 7 would be an example of a 
ropelength-critical knot which is not smooth. 

The only links so far proved to be tight are certain examples from [ |CKS01 |. 
In these links, each component is planar, is piecewise circular (or straight), but fails 
to be unless it is a circle. No explicit tight knots are known. It seems unlikely 
that any of them would be more than piecewise smooth. But our conjectured local- 
minimum trefoil gives the clearest example yet where a discontinuity in curvature 
should happen. 
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Figure 2. Numerically computed curvature (left) and torsion (right) 

functions of the curve 74096, plotted against arclength, going slightly 
more than halfway around the curve. The high frequencies should prob- 
ably be viewed as numerical noise. There is presumably a similar (3, 2)- 
torus knot 7 which is a local minimum for ropelength. The curvature 
of 7 seems to jump abruptly from to i , and then smoothly increase 
further. The torsion is close to zero except where the curvature is high. 
Note that for the finite value p = 4096 shown, the arcs which would be 
straight in 7 are not exactly straight. 
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